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Abstract 

Several studies on real complex networks from different fields as biology, econ- 
omy, or sociology have shown that the degree of nodes (number of edges connected 
to each node) follows a scale- free power-law distribution like P{k) ~ A:"''', where 
P{k) denotes the frequency of the nodes that are connected to k other nodes. Here 
we have carried out a study on scale-free networks, where a line graph transfor- 
mation (i.e., edges in an initial network are transformed into nodes) is applied to 
a power-law distribution. Our results indicate that a power-law distribution as 
P{k) ~ A;"'*'^^ is found for the transformed network together with a peak for low- 
degree nodes. In the present work we show a parametrization of this behaviour 
and discuss its application to real networks as metabolic networks, protein-protein 
interaction network and World Wide Web. 



PACS numbers: 89.75.-k, 05.65.-Fb 

1 Introduction 

Commonly, networks of complex systems were described with the classical random graph 
theory of Erdos and Renyi plEj- In the random graph theory, the degree distribution P{k) 
(probability that a randomly selected node has exactly k edges) peaks strongly around K 
= (k), where () denotes the average. However, recent experimental studies indicated that 
the random graph model could not explain the topological properties of real networks. 
For many real networks the degree distribution was found to follow a scale-free power-law 
distribution like P{k) ~ k^'^ with an exponent 7 between two to four. Some examples 
of these networks are World Wide Web (7 = 2.1), power grid (7 = 4) and film actors 

(7 = 2.3) [Sj. The analyses also showed that some physical networks (neural network 
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of worm C.elegans, film actors, power grid) had common features as small path length 
Lactuai ~ Lrandom but high clustering degree C actual » C random Bi, revealing a different 
topology than the classical random approach. The systems with these properties were 
called " small- world" networks El El • 

Concerning the biological networks, several analyses on metabolic networks corre- 
sponding to 43 different organisms |^ and protein-protein interaction network ^ were 
done. These works revealed the same scale-free topology found in non-biological networks. 

This wealth of data stimulated to develop theoretical approaches to reproduce such 
kind of scale-free topology. One of the most successful models was proposed by Barahdsi- 
Albert jHl which introduced a mean-field method to simulate the growth dynamics 
of individual nodes in a continuum theory framework. The Barabdsi- Albert (BA) model 
is based on two main mechanisms; (1) Growth: the network starts with a small number 
of nodes (mo), and a new node is added at every time step with (m < mo) edges that 
link the new node to (m) different nodes. (2) Preferential attachment: the BA model 
assumes that the probability H that a new node will be connected to node i depends 
on the connectivity ki of that node Y\{ki) = ki/J^jkj. Therefore, after t time steps, the 
network is populated with N = t + niQ nodes and mt edges 0| CHI- Even though recent 
extensions of this model, with rewiring edges , adding a fitness-dependent dynamic 
exponent |T21, and with aging features OClj have provided a more accurate description 
of the network evolution, by generating a large available spectrum of scaling exponent 
or cut-offs in the connectivity distribution ^2j, we will use the original BA model for 
generating our synthetic network. 

The aim of our work is to introduce the concept of the line graph transformation 
and to study the topology of the scale-free networks after the line graph transformation 
is done. As the line graph transformation is closely related to the representation of 
metabolic networks, similarities and differences between the line graph transformation and 
real metabolic networks are also discussed. In addition, we will illustrate our results with 
examples from several real networks as World Wide Web and protein-protein interaction 
networks. 
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It is also worth noting that the hne graph transformation has recently been applied 
with success on the protein interaction network with the aim to detect functional 
modules. In that work, the edges (interactions) between two proteins become the nodes 
of the transformed network (interaction network). By means of the line graph transfor- 
mation, the interaction network has its structure level more increased than that from the 
protein network (i.e., higher clustering coefficient), and by using the TribeMCL algorithm 
PT)] they are able to detect clusters in the more highly clustered interaction network. 
These clusters are transformed back to the initial protein-protein network to identify 
which proteins conform functional clusters. 

The paper is organized as follows. In Sec. II we describe the theoretical model that 
we use and we explain in detail the mathematical methods. In Sec. Ill we present the ex- 
perimental data of several real networks and we compare with our theoretical predictions. 
The final section summarizes our work. 

2 Theoretical Models 
2.1 Line graph transformation 

Given an undirected graph G, defined by a set of nodes V{G) and a set of edges E{G), we 
associate another graph L{G), called the line graph of G, in which V{L{G)) = E{G), and 
where two vertices are adjacent if and only if they have a common endpoint in G (i.e., 
E{L{G)) = {{{u,v),{v,w)}\{u,v) e E{G),{v,w) G E{G)}). This construction of graph 
L{G) from the initial graph G is called line graph transformation |17j . 

In Fig. 1(a), we consider a graph G and we apply the line graph transformation over 
this graph G. The result of this transformation is the line graph of G, L{G). We see 
that, under the line graph transformation, the nodes of L{G) are the edges of G, with 
two nodes of L{G) adjacent whenever the corresponding edges of G are. If x = {u,v) is 
an edge of G, then the degree of x (deg (x)) (i.e., number of edges connected to the node 
x) in L{G) can be written as: deg (u) + deg (v) - 2. The example shown in Fig. 1(a) 
can correspond to a synthetic network generated by the BA model, where the network is 
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only composed by one kind of nodes. In this case, the edges of the initial network E{G) 
do not correspond to any physical entity. The line graph transformation becomes more 
meaningful when those edges E{G) have physical meaning. 

This concept of line graph transformation has similarities to the representation of 
metabolic networks as we can see in Fig. 1(b). Fig. 1(b) (left) shows a real subgraph from 
the Lysine Biosynthesis metabolic network in the ordinary representation of biochemical 
pathways. We see that the nodes are chemical compounds and edges are represented by 
chemical reactions ^. In this example, nodes and edges have physical meaning. There 
are five chemical compounds and four chemical reactions in that representation. /^From 
this representation we can construct the graph G which only contains compounds as 
nodes (chemical compound network). Two nodes are linked by an edge when they occur 
(either as substrates or products) in the same chemical reaction. After applying the line 
graph transformation on the graph G, we obtain the graph L{G) which only contains 
reactions as nodes (reaction network). The graph G (compound network) will have a 
degree distribution Pcik), and the graph L{G) (reaction network), obtained through 
the fine graph transformation, will allow us to study its degree distribution Pl(g){.^) (or 
other topological observable). In particular, we focus on the correlation between both 
degree distributions Pcik) and PL{G){k), through the line graph transformation. In these 
figures, we see clearly the close relationship between the line graph transformation and 
real metabolic networks illustrating the motivation of our analysis. In this example, the 
line graph transformation becomes more meaningful. 

However, we must notice that in some cases the line graph transformation could give 
rise to spurious nodes. To be concrete, we will use the terminology of chemical com- 
pounds/reactions, but this explanation can be extended to any other similar graph struc- 

^It is important to comment the differences between enzymes and chemical reactions. We have con- 
sidered the edges of metaboHc networks as chemical reactions by following the notation from the KEGG 
database. The chemical compounds are also called substrates. Each substrate can be represented as a 
node of the graph, linked by chemical reactions. The products of these chemical reactions appear as 
other nodes (substrates). The enzymes are chemical entities which catalyze the chemical reactions. In 
Fig. 1 we show the Enzyme Commission |18| numbers [EC a.h.c.d\ for several enzymes as an example. 
One enzyme can catalyze more than one chemical reaction. In our work we only focus on the chemical 
reactions, and they are represented by the edges in the metabolic network. 
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ture. These spurious nodes appear when the reactions have more than one product (or 
substrate) as chemical compounds. When these cases appear, the compound graph G 
generated (as explained later with Fig. 6) will have extra edges, and consequently the 
line graph transformation will generate reactions as nodes which do not exist in the real 
metabolic process. We have analysed this issue and a detailed explanation for the case of 
metabolic networks can be found in Sec. III. 




L-2-amino- UDP-N-acetylmuramoyl- L-2-amino- UDP-N-acetylmuramoyl- 
6-oxopimelate L-alanyl-D-glutamyl- 6-oxopimelate L-alanyl-D-glutamyl- 

meso-2,6-diaminopimelate meso-2,6-diaminopimelate 



(b) 

Figure 1: (a) Graph G is an initial network. L{G) is the corresponding line graph network, 
(b) A representation of one subgraph from the Lysine Biosynthesis metabolic network. The 
graph G only contains compounds as nodes (chemical compound network). Two compounds 
are linked by an edge when they occur in the same chemical reaction. Applying the line graph 
transformation we generate the line graph L(G), which is the chemical reaction graph embedded 
in the metabolic network, where the nodes are the reactions. This example illustrates the close 
relationship between the line graph transformation and the metabolic network. 

The main feature underlying the line graph transformation can be summarized as 
follows: We assume that the initial network has the scale-free topology as P{k) ^ k~'^. 
As the degree of each transformed node (i.e., an edge in G) will be roughly around k, 
the distribution of the line graph L{G) should be k ■ k~"' = k~'^^^ with degree around 
k. Therefore, we can conclude that if we have a graph G with a probability distribution 
following a power-law as k''' , then L{G) will follow a power-law as k~"'~^^. We have de- 
veloped two models to reproduce the behavior of the line graph transformation over the 
scale-free network. In the first one we solve the discrete equation for the degree distri- 
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bution of the line graph and in the second one we use the inverted beta distribution. A 
detailed mathematical explanation on these models can be found in the next subsections. 



1 1 1 1 — I — I — r~r"| 1 1 1 1 — I — I — r~r"| 1 1 1 1 — i — i — tte 
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Figure 2: Eq. (3) is numerically computed with the BA model by fixing degree j = jo (jo = 
10, 20, 50). 30 trials of BA networks of 10000 nodes are computed for each jo in order to have 
enough statistics. The results reproduce Eq. (4) with 7' = —2 (continuous lines). 



2.2 Discrete equation for the degree distribution of a trans- 
formed graph 

We asumme the following: A) degree distribution of an original scale-free network follows 
\{v\deg{v) = oc if {d > mo), and \{v\deg{v) = d}\ =0 otherwise. B) edges are 
randomly generated under (A). Precisely, deg{u) is independent of deg{v) for each edge 
(m, f). Hence, if deg{u) is 1 + di and deg{v) is 1 + ^2, the transformed node correspond- 
ing to edge ("U, f) has degree di + d2. The following relations hold from {A) and {B): 
Proh [degiv) = 1 + rfj] oc (1 + di)'"' with i=l,2, for a randomly generated edge {u, v), 
where di + 1, d2 + 1 > rriQ. Therefore, degree distribution of a transformed node would 
be given by: 
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Prob[deg{u,v) = d](x ^ (1 + (^1)-^+^ ■ (1 + ^2)"^+^ (1) 

d=di+d2 

This equation can be expressed as: 

do -mo 

Prob [deg{u, v) = d] <x ^ [x^-^ ■ {do - x^'^'] (2) 

x=mo 

with do = d + 2 and x = 1 + di. We are also able to sum in a discrete way in terms of 
polygamma functions (see the next subsection). 

Though we have assumed that deg{u) is independent from deg{v), this assumption 
can be weakened to the condition that the original network shows no assortative mixing 
(neutral network). By assortative (disassortative) mixing in networks we understand the 
preference for nodes with high degree to connect to other high (low) degree nodes fK)\ E] . 
We could choose randomly one edge and we consider the node reached by following that 
edge. Following Newman the degree distribution for the node at the end of a randomly 
chosen edge will be: {k + 1) ■ P{k +1), if we only consider the number of edges leaving the 
node (i.e., not taking into account the node we arrived along). Therefore, the probability 
distribution ejk of all edges {u,v) that link together nodes with degree j + k (sum of the 
degrees of the nodes at the ends of (m, v) edge) would be approximated as: 

e,fc oc (j + 1) ■ (j + l)-^k + l)-ik + (3) 

In order to validate Eq. (3) with the BA model, we compute ejk by using several fixed 
values of j = jo- By fixing j, Eq. (3) can be written as: 

e,„fcCxC(jo)-(A: + l)-^+\ (4) 

where C{jo) is a constant, and as 7 = 3 for the BA model, 7' = —7 + 1 should be —2. We 
show the results of our numerical computation in Fig. 2. We see three power-law with 
exponent 7' = —2 for three different values of jo- These results indicate that Eq. (3) is a 
valid expression for the BA model. It is worth noticing that ejk = Ckj- 
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Figure 3: Analysis with synthetic network, (a) Circles (blue): degree distribution generated by 
using the Barahasi- Albert model with mo=m=5 (in what follows, mQ=m). The synthetic 

network follows the power-law distribution P{k) ~ k~'^ with exponent 7 = 2.9 it 0.1 (continuous 
line). Diamonds (red): the distribution of the transformed scale- free network also follows a 
power-law with exponent 7= 1.9 it 0.1 (dashed line). We took average of 10 BA networks of 
10000 nodes. (b) Inverted beta distribution (continuous line) with parameters /3 = 17, a = 1 and 
a = rriQ = 5 and polygamma distribution for 7 = 3 (dash-dotted line). In the inset we show 
the transformed network data together the polygamma distribution in linear-linear scale. It is 
interesting to note that the peak found in the transformed network does not exist in the original 
network. This feature can also be predicted by our theoretical model. 

In order to compare with Eq. (1), we sum for all the possible degrees of the two 
vertices at either end of edges {u,v) which can generate transformed nodes with degree d: 



Prob[deg{u,v) = d]= ^ e.^ oc Yl (j + 1)"^+' • (A: + l)-^+\ (5) 

d—j'\-k d—j-\-k 
j + l>m0 J-|-l>mo 
fc + l>m0 fc + l>mQ 

and Eq. (1) holds. 

Although in [201 it is suggested that the BA model is not assortative (neutral scale-free 
network) (i.e., Pearson correlation function gives a zero result), Redner and Krapivsky [221 
pointed out that obvious degree correlations exist in BA model (older nodes tend to be 
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connected) and in similar growth models. However, as we have seen in Fig. 2, Eq. (3) is 
a good approximation for the BA model. 

2.3 Analysis in terms of Poly gamma functions 

The digamma function is defined as: ^{z) — ^lnr{z) — where r{z) is an extension 
of the factorial to complex and real numbers and it is analytic everywhere except at 
z — 0,-1,— 2,— 3.... The nth derivate of ^ is named the Poly gamma function, denoted 
ipn{z)- The notation ipo{z) is frequently used for the digamma function itself as ipniz) — 
-^Y^ = ^^o(-2)- Other useful identities are 

dz + 

and 

V'o(x) = -7 + Et (7) 

fc=i 

where 7 is the Euler-Mascheroni constant and x is a positive integer value. 

By decomposing Eq. (2) in partial fractions, and by using the polygamma relationships 
shown above, we can find the expressions for Eq. (2) evaluated with 7 = 2 and 7 = 3. 
We detail the calculation for 7 = 2. The summation of Eq. (2) can be decomposed in 
two terms: 

do-mo ^ S 

where A = B = l/d^. 

The first term of the sum gives: 

do-mo 111 1 1 

E - = — + 7 + - + 3 7 + 3 9 

^t;^^ X mo mo + 1 do - mo + 1 do - mo 

and by using Eq. (7), we can write it as: ipo{x + 1) — '^o(^o)) with x = do — mo- 
Analogously, the second term gives the same contribution. Therefore for 7 = 2, Eq. (2) 
reads as: 
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Prob [deg{u, v) = d\^=2 oc —[ipoix + 1) - V^o("^o)] 



(10) 



For 7 = 3 we have 



Prob [deg{u, v) = d\'y=3 oc 



(11) 



with 



ipo = 2[iJo{x + 1) - i^o{mo)] 



(12) 



and 



V^i =2[^i(mo) -V^i(a; + 1)], 



(13) 



where do = d + 2 and x = do — mo are integer values. 

It is interesting to note that Prob[deg{u,v) = d]^=2 ~ d~^ and Prob[deg{u,v) = 
d].y=3 ^ for large d, which matches the distribution of k~'^^^ as we can see in Fig. 



2.4 Inverted beta distribution 

One drawback of the previous approach is that the exponent 7 is considered as an integer 
number which constrains its range of quantitative applicability. Hence, we have looked 
for a continuous function in terms of the exponent 7. In that sense, we have found 
that the inverted beta distribution B{y) oc {y — a)^^^/{l + y — a)"^^, which is obtained 
doing the transformation Y = over the beta distribution followed by a translation 
y ^ y — a, fits well our requirements and it was used successfully to reproduce the data. 
In addition, it is interesting to note that the distribution shows a power-law tail for large 
y as: B{y) y"""^ 

2.5 Theoretical results 

Once the theoretical approaches have been introduced, we generate a synthetic scale-free 
network with exponent 7 = 2.9 using the BA model [HI CD] and we study the behavior of 



3(b). 
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the line graph transformation over that scale-free network. We compare the transformed 
network with the results from the theoretical functions introduced in the previous sections. 
In Fig. 3(a) we see that the transformed network (diamonds) follows a power-law with 
exponent 7 = 1.9. It is interesting to note that the exponent of this scale-free is exactly 
reduced by one unit (from the original scale-free network to the transformed network) as 
it was predicted in the previous sections. As a second result, a peak was found for low 
degree nodes in the line graph transformed network, indicating that the power-law is like 
a tail (or asymptotic behaviour) of a more general distribution. In Fig. 3(b) we plot the 
inverted beta distribution and the polygamma function to compare with the transformed 
network. We see that the curves reproduce well the peak of the data for low degree k 
and also shows a power-law tail for higher degree k. Both agreements give us confidence 
about the fairness of both approaches used to study the data. 

3 Experimental data 

There are several examples, in biological and non-biological networks, which appear to 
support our theoretical analysis. We have tested our results in the World Wide Web 
network with a size of 325729 nodes representing web pages being connected by links 
each other. The dataset was obtained from the website of Notre-Dame Research Group 
[SB]- We must notice that we have considered the links as undirected edges in order to 
compare with our theoretical approach. In that sense the value of 7 obtained here could be 
considered as an average of the 7j„, •jout [IH]- We have also analysed the protein-protein 
interaction network for the yeast S. Cerevisae which contains around 1870 proteins as 
nodes linked by 2240 bidirectional interactions [HII2S1- In Fig- 4 and Fig. 5 we show the 
data for the WWW network and the protein-protein interaction network respectively. We 
see that both networks (circles) are following the power-law P{k) ~ k~'^. In the same 
figures, we present our results for the corresponding transformed network (diamonds). In 
both cases, we see that the scale-free topology is preserved and the exponent 7 is decreased 
by almost one unit as we expected. 
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Figure 4: Circles (blue): degree distribution for the VKM^VK network generated with data from 
j23j . The VKM^VK network follows a power-law with exponent 7 = 2.27. Diamonds (red): the 
distribution of the transformed network also shows a power-law tail with exponent 7 = 1.32. 
We show the inverted beta distribution evaluated with o = mo = 1 (dashed line), a = ttiq = 
(dash-dotted line), and the beta parameters are /3 = 16, a = 0.4. The correlation coefficient is 

Following with our analysis of real networks, we analyse now the metabolic network 
from the KEGG database fI^ which contained around 10400 chemical compounds and 
4100 chemical reactions. The KEGG database is one of the best sites for biochemistry, 
metabolism, and molecular biology information. As we showed in Fig. 1, the line graph 
transformation is closely related to the representation of metabolic networks, however we 
must point out some differences which we have found in our study. 

As it is depicted in Fig. 6, in the metabolic networks exist some cases where not all the 
edges associated to the substrate network should be transformed by the line graph trans- 
formation. The main issue is that for each chemical reaction with more than one product 
(or substrate), we must only transform the same number of edges from the substrate graph 
as the number of reactions in the real metabolic process. In the case shown in Fig. 6, the 
resulting network (reaction network), after an ordinary line graph transformation, would 
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Figure 5: Circles (blue): degree distribution for the protein-protein interaction network S. 
Cerevisae generated with data from The protein-protein interaction network follows a 

power-law with exponent 7 = 2.20. Diamonds (red): the distribution of the transformed network 
also shows a power-law tail with exponent 7 = 1.32. We show the inverted beta distribution 
evaluated with a = uiq = 1 (dashed line), a = niQ = (dash-dotted line), and the beta 
parameters are /? = 6, q = 0.4. The correlation coefficient is r^. 

generate up to nine nodes (reactions), but only two exist in the real metabolic process. 
These spurious nodes only appear when one (or some) reaction(s) in the network have 
more than one product (or substrate). These cases should be transformed as it is shown in 
Fig. 6 (i.e., by creating only as many nodes in the transformed network as reactions in the 
real metabolic process). This procedure is called physical line graph transformation. The 
graph generated by the physical line graph transformation is denoted by L^eaiiG)- In a 
network where the are some spurious cases as exposed above, we proceed by applying the 
physical line transformation for these ones, and the ordinary line graph transformation 
for the rest of the cases. The resulting distribution could be different if we compare it 
with the distribution generated only by the line graph transformation P{k) ~ k"^'^^ . The 
distortion could be larger if there are many of these spurious cases in a network. 

In Fig. 7 we show the degree distribution of the chemical compounds (circles) in the 
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Figure 6: Sketch of metabolic networks. We show two reactions {Rl, R2), with only one 
common substrate (product) C4. The substrate graph G (compounds) is shown with dark blue 
circles. The reaction graph Lreai{G) (reactions) is shown with light red circles. To obtain this 
graph L^eaiiG) we apply a physical line graph transformation, i.e., we do not transform all the 
edges from the initial graph and we only generate the same number of nodes in the transformed 
network as the number of reactions in the real metabolic process. An ordinary line graph 
transformation would generate nine reactions as nodes in the line graph L{G), when only two 
reactions exist in the real network. 

metabolic network from the KEGG database. In the same figure, we have plotted the 
distribution (diamonds) for the transformed network which corresponds to the reaction 
network. We see that both distributions follow a power-law and the difference between 
their exponents is one. However, that transformed network L{G) would not have a full 
physical meaning because it could have extra reactions as nodes which do not exist in 
the real metabolic process, as we explained in Fig. 6. In Fig. 8 we plot the transformed 
network after applying the line graph transformation and the physical line graph trans- 
formation for the spurious cases. We see that both graphs are scale-free networks, but in 
this case, the difference between the exponents is smaller than one, due to the distortion 
mentioned above. 

Experimental results shown in our paper suggest that exponents of the transformed 
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scale-free network may change due to the effect of the spurious cases. However, the 
experiments also indicate that the scale-free topology is preserved. In particular, if the 
average number of substrates and products per reaction is small, the spurious cases may 
not strongly affect, and the scale-free power-law distribution would be preserved. 




k [Degree] 



Figure 7: Circles (blue): Degree distribution for the metabolic networks of KEGG database 
j24j . The data included the following 9 metabolic pathways: Carbohydrate Metabolism, Energy 
Metabolism, Lipid Metabolism, Nucleotide Metabolism, Amino Acid Metabolism, Metabolism of 
Other Amino Acids, Metabolism of Complex Lipids, Metabolism of Complex Carbohydrates and 
Metabolism of Cof actors and Vitamins. Diamonds (red): The metabolic network follows a power- 
law with exponent 7 = 2.24 (continuous line). Diamonds (red): the transformed distribution 
of the scale-free network shows a power-law tail with exponent 7 = 1.21 (dashed line). With 
dash-dotted line, we show the inverted beta distribution with parameters /3 = 2.5, a = 0.25 
evaluated with a = rriQ = 1. The correlation coefficient is r^. 

It is interesting to quote a previous work [23], which analysed the distribution of 
metabolite connectivities in both substrate and reaction graphs. Although in that paper 
they do not mention about the line graph transformation technique and, consequently, 
they do not discuss about the reason of the modification of exponent 7, they notice that 
the degree distribution in the reaction graph does not follow a simple power-law and it 
appears to be governed by two quantitatively different regimes. 
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Figure 8: Circles (blue): Same as Fig. 5. Triangles (green): the distribution of the transformed 
network after applying the line graph transformation and the physical line graph transformation 
for the spurious cases. It shows a power-tail with exponent 7 = 1.70 (dashed line). With dash- 
dotted line, we show the inverted beta distribution with parameters P = 2.7, a = 0.7 evaluated 
with a = mo = 1- The correlation coefficient is r^. 

4 Conclusions 

We have reported on the two complementary representations of a scale-free network using 
the line graph transformation. This transformation is useful when it is applied on networks 
where the edges have physical meaning. In particular, we have illustrated that the line 
graph transformation is closely related to the representation of metabolic network. In this 
network by using the line graph transformation, the reaction network can be generated 
and it allows us to study its degree distribution P{k) and other topological observables 
as clustering degree, for example. 

The two goals of the present work can be summarized as follows: We have described 
the real networks as two complementary representations of a scale-free network, where 
the second one emerges when the line graph transformation is done over the first one. 
Our second goal is that we have found that the exponent 7 is always one unit less than 
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the initial exponent 7 coming from the original scale-free network. We have seen that 
the difference of one unit in the exponents of the degree distributions can be found in 
biological and non-biological networks. We have also carried out a theoretical study of 
the general distribution underlying the line graph transformation, being successfully to 
reproduce the tail of the power-law and the peak found for low-degree nodes. However, we 
should bear in mind that in some particular cases (e.g., metabolic networks) an ordinary 
line graph transformation applied over a real network could generate extra nodes in the 
transformed network without a real correspondence. These spurious cases make that the 
difference between the exponents 7 could differ from one unit. 

We also took advantage of the fact that the BA model is not assortative 120] to 
compare our results for the discrete equation of a transformed graph. As a future work, 
an extension of this analysis to assortative networks should be done. 

This study is an interesting step forward to understanding large complex networks 
from this complementary scale-free perspective. 
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